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SUMMARY 


The  problem  of  two-dimensional  flow  behind  a  curved  stationary 
shock  wave  is  considered  analytically.  The  method  assvm.es  a  given 
shock-wave  shape,  which  automatically  determines  certain  initial  con¬ 
ditions  on  the  flow  variables;  and  the  flow  pattern,  including  any 
body  shape,  follows  from  the  initial  conditions.  Approximate  analytic 
expressions  are  found  for  the  stream  function  in  the  subsonic  region 
following  the  shock  and,  after  the  stream  function  is  obtained,  the 
flow  density  is  determined  by  Bernoulli's  equation  which  connects  the 
density  with  the  derivatives  of  the  stream  function.  The  final  solution 
can  then  be  determined  from  the  velocity  field  thus  obtained. 


INTRODUCTION 


The  problem  of  compressible  flow  after  a  curved  shock  wave  has  been 
investigated  by  various  authors  (references  1  to  4).  In  the  analytical 
treatments,  the  main  interest  has  so  far  been  concerned  with  the  local 
properties  of  the  flow,  such  as  the  relations  between  the  gradients  of 
various  physical  and  geometrical  quantities  along  the  shock  wave  and 
those  along  the  body. 

The  purpose  of  the  present  work  is  an  attempt  to  treat  analytically 
the  two-dimensional  problem  after  the  shock  in  the  large.  The  method 
assumes  a  given  shock-wave  shape,  which  automatically  determines  certain 
initial  conditions  on  the  flow  variables.  It  is  therefore  a  Cauchy 
problem  and  the  flow  pattern,  including  any  body  profile,  follows  from 
the  assumed  initial  conditions.  Approximate  analytic  expressions  for 
the  stream  function  are  found  for  the  compressible  flow  in  the  whole 
subsonic  region.  The  stream  function  thus  found  satisfies  the  exact 
shock  conditions  but  it  satisfies  the  differential  equation  only  approxi¬ 
mately.  After  the  stream  function  is  obtained,  the  flow  density  is 
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determined  by  the  exact  Bernoulli  equation  which  connects  the  gas 
density  with  the  derivatives  of  the  stream  function.  Then  the  sonic 
line  can  be  determined  from  the  velocity  field  thus  obtained. 

The  author  wishes  to  take  this  opportunity  to  thank  Professor 
C.  C.  Lin  for  his  valuable  criticisms  and  discussions. 

This  work  was  conducted  at  the  Massachusetts  Institute  of  Technology 
under  the  sponsorship  and  with  the  financial  assistance  of  the  National 
Advisory  Committee  for  Aeronautics. 


SYMBOLS 


7 

Pi 

Pi 

u 


Cl 

F(PiUy) 

(xo>0) 


P.Q 


ratio  of  specific  heats  of  fluids 
free- stream  pressure 
free- stream  density 

free-stream  velocity  magnitude,  assumed  to  be  parallel 
to  x-axis 

free-stream  velocity  of  sound 
function  determining  shock  shape  (l) 

shock  nose  when  F(0)  =0 
stagnation  density  after  shock  wave 
stagnation  velocity  of  sound  after  shock  wave 

where  \|r  is  stream  function 
complex  variable  in  holograph  plane 

((♦*  -  *  ‘’o'o)) 

parameteric  representation  of  a  three-dimensional  curve 
with  parameter  r 

complex  variables 
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equation  for  shock 


Ml 

tan  T 
P 
P 

G(ik) 

F(t) 

Z,m 
u,v 

q  =  \|u2  +  v2 

c  local  velocity  of  sound  after  shock  wave 

maximum  velocity  magnitude  of  flow 

s,n  tangential  and  normal  directions  at  point  on  shock  wave 

q^,qj^  tangential  and  normal  components  of  velocity  of  flow 

immediately  after  shock  wave 

K2_^K2  constants  defined  hy  equations  (l8) 

F(y)  analytic  function  in  real  variable  y 

0(x,y)  auxiliary  function  defined  by  equations  (22) 

h{Q>)  arbitrary  analytic  function  of  Q 

Re  real  part  of  a  complex  quantity 


free- stream  Mach  number  ((TPi/Pij-Z-U 

stream  function  of  flow  after  shock  wave 
the  shock) 

slope  of  shock  wave  at  point  on  it 
pressure  after  shock  wave 
density  after  shock  wave 
entropy  function 
function  determining  shock  shape 
constants  defined  by  equations  (8) 
velocity  components  of  flow  after  shock  wave 


h 
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k  =  PiU/p^Cq 


z  =  X  +  yi 

fi(t) 

t 


Im 


i 


n(P).t(Q) 


K(r) 

n(y) 

cr  =  y  +  io) 
0) 


analytic  fTinctlon  for  Weierstrass’  parametric  represen¬ 
tation  of  a  minimal  surface 

complex  variable  occurring  in  Weierstrass^  formula  for 
minimal  surfaces 

imaginary  part  of  a  complex  quantity 

vector  function  defined  by  equation  (29) 

vector  with  components  x,  y,  and  w 

complex-valued  vector  functions  of  variables  P  and  Q, 
respectively;  components  denoted  by  subscripts  1,  2, 
and  3 

unit,  vector  normal  to  tangential  direction  of 

vector  function  defined  ty  eq\iation  (33) 

complex  variable 

imaginary  part  of  complex  variable  a 


FORMULATION  OF  PROBLEM 


The  present  report  is  concerned  with  a  stationary  shock  wave  pro¬ 
duced  by  a  two-dimensional  body  placed  in  a  tinlform  supersonic  stream. 
The  fluid  is  assumed  to  be  an  ideal  gas  and  the  process  on  any  material 
particle,  to  be  adiabatic.  Thus,  the  entropy  is  constant  along  a 
streamline  while  it  may  vary  from  one  streamline  to  the  next. 


Let  the  coordinate  axes  be  so  chosen  that  the  incoming  free  stream 
is  along  the  x-axis  and  that  the  shape  of  the  shock  wave  is  represented 
by 


'  [5(pi'Jy)/pioJ 


+  X. 


X 


(1) 
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where  U  is  the  free- stream  velocity  and  p-j_  is  the  gas  density  in 
the  free  stream.  As  the  shock  wave  must  asymptotically  become  the 

Mach  line  at  infinity^  the  function  F*  (p^^y)  — >  -  1  as 

y  — >  ±09,  where  is  the  free- stream  Mach  number. 

Entropy  of  flow  after  shock  waves.-  The  stream  function  ilr(x,y) 
may  be  regarded  as  the  amount  of  the  flow  across  a  curve  joining  the 
point  (x,y)  to  the  point  (xq,o).  It  must  therefore  be  continuous  across 

the  shock  wave  because  of  the  continuity  of  the  mass.  Thus,  on  the 
shock  wave, 


t(x,y)  =  p^Uy 


(2) 


The  slope  of  the  shock  wave  is  therefore 


tan  T  =  ^  =  l/F'(\lf) 


Since  the  entropy  is  constant  along  a  streamline,  the  pres sure- density 
relation  for  the  flow  after  the  shock  wave  can  be  written  in  the  form 


p/p’'  =  G(\|f)p^yp^’' 


(4) 


where  y  is  the  ratio  of  the  specific  heats  of  the  fluid,  -p^  and  p 
are,  respectively,  the  pressures  before  and  after  the  shock,  and 
and  p  are  corresponding  densities  of  the  fluid.  This  entropy  function 
G(t)  is  now  to  be  expressed  in  terms  of  the  function  F(^lf)  which  deter¬ 
mines  the  shock  shape.  From  the  shock  conditions,  it  can  be  immediately 
shown  that 


p/pi  = 


27Mj^  sln^T  +  (l  -  7) 
7+1 


p/pi  = 


(7  +  i)Mi^  sin^T 
2  +  (7  -  sin^T 


(6) 
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it  follows  that 


I 


.2ryy 


^  g' ii) 

27  G(i) 


2  2 
P  '3m 


1  + 


(12) 


where 


+  V' 


(13) 


c2  =  7p/p  (l4) 


and  qj^  is  the  maximum  velocity  of  the  flow.  This  eqmtion  is  here¬ 
after  referred  to  as  the  fundamental  equation  for  compressible  flow 
after  the  curved  stationary  shock  wave  given  hy  equation  (l).  In  it^ 
the  quantities  u,  v,  c,  and  p  are  known  functions  of  \[f^, 

and  determined  hy  equations  (9)^  and  (ll). 

It  should  he  noted  that  equation  (l2)  is  equivalent  to  Crocco’  s 
equation  (reference  5)  for  vortex  motion.  In  fact,  one  is  immediately 
derived  from  the  other  hy  the  relation  between  the  stream  function 
defined  above  and  that  defined  hy  Crocco. 


Shock  conditions  expressed  in  terms  of  derivatives  of  stream 
function. -  Let  s  and  n  he,  respectively,  the  tangential  and  the 
norma l^ire c t ions  at  a  point  on  the  shock  wave.  Then 


cos 


Ss 


T  + 


Sin  T 


Sn 

Sn 


sin  T 


cos  T 


(15) 


From  continuity  of  the  flow  across  the  shock  it  follows  that 


g  =  pq^  =  PiU  sin  T 

^  TT 

^  =  -pq^  =  -pU  cos  T 
on 


(16) 
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where  and  are,  respectively,  the  tangential  and  the  normal 

components  of  the  velocity  of  the  flow  immediately  after  the  shock 
wave.  Thus,  by  virtue  of  equations  (15)  and  (6),  the  shock  conditions 
in  terms  of  the  derivatives  of  the  stream  function  are  transformed  into 
the  following  form: 


Mathematical  formulation  of  problem.-  The  problem  of  finding  the 
two-dimensional-flow  pattern  after  a  stationary  shock  wave  given  by 
equation  (l)  is  then  to  find  the  solution  of  the  fundamental  equa¬ 
tion  (12)  with  the  conditions  (17)  on  the  shock  wave.  It  is  therefore 
a  Cauchy  initial-value  problem  if  the  shock  wave  is  given. 

However,  it  should  be  noted  that  the  stream  function  >K(x,y)  thus 
obtained  is,  in  general,  many-valued  if  it  is  extended  continuously  (in 
physical  variables)  into  the  whole  region  after  the  shock  wave.  In  the 
actual  case,  it  is  expected  that  other  shock  waves  are  formed  near  the 
tail  of  the  body.  Therefore  the  solution,  once  obtained,  is  limited 
by  such  shocks.  The  flow  farther  downstream  should  be  treated  separately. 
In  fact,  the  flow  after  the  shock  wave  can  be  treated  separately  for  the 
subsonic  region  and  the  supersonic  region. 

Determination  of  physical  variables.-  After  t(x,y)  is  found,  the 
density  distribution  may  be  obtained  from  equations  (4)  and  (ll).  The 
velocity  components  may  then  be  gotten  from  equation  (9)*  It  should  be 
remarked  that  this  determination  is  two-valued:  One  value  of  the 
density  corresponds  to  supersonic  flow  and  the  other,  to  subsonic  flow. 

The  location  of  the  sonic  line  is  therefore  essential.  When  a  solution 
is  obtained  for  t  (whether  exact  or  approximate)  the  sonic  line  is 
immediately  given  by  the  condition  that  the  two  roots  in  Bernoulli's 
equation  for  the  density  p  are  eqiial.  This  gives 


/  p  y-1  ^  7-1/2  2^  ^  ^  1 

7  +  l\7  +  1  ^  /G(i|r)  G(\|f) 
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Substituting  back  into  Bernoulli' s  eqmtion,  it  follows  that 


where 


(18a) 


(l3b) 


and  c-j^  is  the  free- stream  velocity  of  soimd. 

In  order  that  the  sonic  line  thus  determined  should  pass  through 
the  sonic  point  immediately  after  the  shock  wave,  must 

satisfy  the  exact  shock  condition  on  the  shock  wave. 

Approximations,-  In  the  following  discussions,  restrictions  are  put 
on  the  shock  shape  such  that 

(1)  The  shock  has  its  nose  at  (xq,o)  where  the  shock  is  normal; 
that  is,  f' (0)  =  0. 

(2)  F(y)  is  analytic  in  the  real  variable  y;  moreover,  F"(y)  is 
positive. 

(3)  The  shock  wave  tends  to  a  Mach  line  asymptotically. 

From  the  first  two  assumptions  it  follows  that  F' (y)  Is  a  monotonic 
function  of  y  and  yF'(y)  ^  0.  From  the  last  asstimption,  it  is  seen 

readily  that  0  $F'(y)  =  ^M^^  -  ij  ^  . 

The  effect  of  vorticity  generated  by  the  shock  can  be  better  under¬ 
stood  by  examining  the  term  in  the  right-hand  side  of  the  fundamental 
eqioation  (12)  for  From  equations  (3),  the  term  on  the  right-hand 

side  of  equation  (12)  is  a  product  of  F"(A|f)p^qj^^  and  a  nondimens ional 
quantity.  By  a  straightforward  calculation  this  nondimens ional  factor  is 


F’(i) 

fT 

r 

1+  (F’  )2 

m2_1 

L_ 

V 

r 

7 

1  27  _ 

2  1+  (f'  )^ 

m/  +  2  - — 4£_J_ 

7  “  1  _ 
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A  rough  estimation  shows  that  it  is  bounded  by 


2  -lllllLl— _J_ 

(Mi^  _  i)l/2  1  +  (F' 


It  is  therefore  clear  that,  for  the  free- stream  Mach  number  close 

to  1,  a  perturbation  process  based  on  the  parameter 


can  be  performed.  The  term  in  the  right-hand  side  of  equation  (12) 
can  be  neglected  for  the  first  approximation.  In  fact,  for  any  Mach 
n-umber  M^^  in  the  general  case,^  the  nondimensional  factor  is  an 

infinitesimal  of  first  order  for  small  values  of  and  of  second 
order  for  large  values  of  The  assumption  of  neglecting  the  effect 

of  vorticity  is  perhaps  reasonable  for  a  certain  range  of  M2_  greater 
than  unity. 

In  the  following  discussion  the  effect  of  vorticity  will  be  neg¬ 
lected  in  finding  approximate  solutions  for  the  stream  function.  The 
subsonic  region  will  now  be  treated  analytically  by  finding  approximate 
solutions  for  the  stream  function.  The  exact  shock  conditions  will  be 
satisfied  but  the  f-undamental  eqiiation  (12)  will  be  satisfied  only  with 
the  following  approximations: 

(1)  The  effect  of  vorticity  is  neglected 

(2)  Equation  (12)  is  simplified  by  an  approximate  estimation  of  u/c 
and  v/c  as  explicit  functions  ^  and 

This  estimation  follows  the  same  line  as  the  one  first  suggested  by 
Chaplygin  (reference  6).  The  idea  has  been  successfully  applied  by 
many  authors  (references  7  to  10)  for  analytic  treatment  of  subsonic 
flows. 


^An  upper  bound  for  the  nondimensional  factor  in  terms  of  M^^  only 
can  be  determined.  The  following  is  a  trivial  one: 
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It  may  further  be  suggested  that  the  solution  of  the  simplified 
equation  for  satisfying  the  exact  shock  conditions  may  also  serve 
as  a  first  approximation  for  the  supersonic  region  in  the  case  of  a 
thin  body  with  a  blunt  nose.  For  such  a  case,  the  streamlines  are 
likely  to  be  nearly  straight  in  the  supersonic  region.  As  the  approxi¬ 
mate  solution  of  i|r  obtained  from  the  above  considerations  represents 
the  true  geometrical  pattern  of  streamlines  on  the  shock  wave,  far  away 
from  the  body,  and  near  to  it,  it  may  be  expected  to  be  a  reasonable 
approximation  in  the  whole  region.  However,  this  can  be  applied  only 
if  the  two-dimensional  body  is  an  analytic  curve.  In  case  the  body  has 
a  discontinuity,  the  flow  downstream  should  be  treated  separately  from 
the  Mach  line  passing  through  this  point. 

After  having  obtained  the  first  approximation  (covering  both  the 
subsonic  and  the  supersonic  regions),  an  iteration  procedure  may  then 
be  carried  out.  This  is  done  by  evaluating  all  quantities  in  equa¬ 
tion  (12)  except  ^yy  approximate  solution 

and  then  solving  the  linear  equation  for 


AN  ANALYTIC  APPROXIMATE  SOLUTION  IN  THE  LARGE 
FOR  SUBSONIC  REGION 


Solution  of  simplified  differential  eqmtlon.-  The  solution  of 
equation  (12)  satisfying  the  initial  conditions  (17)  on  the  shock  will 
now  be  approximately  obtained  by  finding  the  stream  function  1^(x,y) 
which  satisfies  the  following  simplified  equation: 


\f>o^o) 


+ 


XX 


2„  2 


1  - 


,Po^O 


^yy  =  0 


(19) 


with  the  same  Initial  conditions  (l?),-  Equation  (19)  can  be  written  as 


where 


a 

r  ^  ‘1 

r  1 

dx 

Sy 

f-  -  (e^Po^Co^) 

(20) 


(21) 
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Equation  (20)  implies  the  existence  of  a  function  0  such  that 


]jl  -  (^/po^Cq^) 


(b^/PoV) 


(22) 


It  is  now  readily  seen  that  the  above  development  is  completely  equiv¬ 
alent  to  the  Chaplygin  theory.  However,  the  usTial  physical  interpre¬ 
tations  of  0  and  the  partial  derivatives  of  0  and  ^  are  not 
implied.  The  formal  relations,  however,  remain  ^lnchanged. 

It  is  well-known  that  the  function  0  satisfies  the  eqmtion  for 
the  minimal  svirface  and  that  the  solution  for  0  and  can  be 
expressed  as 


x|f  -  10  =  2pQCQh(fi) 


(23) 


dz  =  dx  +  1  dy 


to 

Q 


n  dh 


(24) 


where  h(fi)  is  an  arbitrary  analytic  function,  and 


Q  = 


2^  2 


- 


Po^o 


(25) 


To  satisfy  the  initial  conditions  (l7),  it  is  fo\ind  convenient  to  apply 
the  theory  for  minimal  surfaces  developed  by  Weierstrass,  Bjorllng, 
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and  Schwarz  (see  reference  ll).  The  final  solution  may  he  given  in  the 
following  differential  form: 


where  Re  denotes  the  real  part  of  a  complex  q\iantlty, 

k  =  PiU/PqCq 

and 

F'  =  F'  (pj^UP) 

and  P  is  a  complex  running  variable.  It  may  he  noted  that  this  solu¬ 
tion  involves  only  the  function  F'(pnUy)  for  the  slope  of  the  shock 
wave.  '  ' 
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Weierstrass^  formula  for  minimal  surfaces^-  Weierstrass^  fomula 
for  minimal  surfaces  (reference  11 )  reads  as  follows: 


where  f^^  is  an  analytic  function  of  the  complex  variatile  t. 

Now 

t  I  dt  =  dh  (l/it)  (28) 

From  equation  (27)  the  following  relation  is  immediately  derived: 
dz  =  dx  +  i  dy  =  it  dh  (l/it)  +  dh  (l/it)/it 


Since  dw  =  Re  itf2”^(t)  dt  =  Re  2i  dh  (H)  =  the  above  rela¬ 

tion  reduces  to  equation  (24)  if  one  sets  =  l/it* 


Shock  conditions  and  Bjorling  problem.-  As  far  as  differential 
equation  (19)  is  concerned,  the  parametric  function  h(n)  in  expres¬ 
sion  (24)  always  generates  an  exact  solution.  There  are  still  shock 
conditions  to  be  satisfied.  These  conditions  (17)  raise  the  problem  of 
finding  the  solution  in  the  large  of  the  Cauchy  initial-value  problem 
for  minimal  surfaces.  Geometrically,  it  is  now  reduced  to  the  problem 
of  finding  the  minimal  surface  when  a  piece  of  space  curve  and  the 
tangent  plane  at  every  point  on  the  curve  are  prescribed.  The  curve  is 

given  in  ?(x,y,w)  space  with  the  following  parametric  representation 
in  vector  form 


y,  0/poCo 


(29) 
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where 


0/poCo  =  0 


F (piUy) 


=  J'  [0x^'(plUy)  +  0^dy 


with  0x  and  0y  determined  by  equations  (22)  and  (l7).  The 
direction  of  the  normal  of  the  tangent  plane  at  a  point 

jF(piUy)/piUj  y,  is  given  by  (0x^0yj-Po'=o)*  This  geometrical 

problem  was  first  suggested  by  Bjorling  in  l844  and  solved  by  Schwarz 
in  1874, 

\ 

Schwarz  fomula.-  Let  ^  he  a  vector  with  components 

and  w.  Let  t](p)  and  ^  (  Q)  he  complex-valued  vector  functions  of 
the  complex  variables  P  and  Q,  respectively.  The  Schwarz  formula  for 
the  Bjorling  problem  with  the  given  curve  |  normal  of 

the  tangent  plane  at  any  point  on  the  curve  (t(r)  such  that  k  =1 
and  K  X  =0,  where  is  the  tangent  vector  of  the  curve,  is 

(reference  ll): 


where 


2?  =t(P)  +  ?(Q) 


l^(P)  =  dtQ(p)  +  i[K{p)  X  d'fo(P)] 
L?(Q)  =  dl^CQ)  -  ipCQ)  X  df^CQ)] 


From  the  last  two  relations  it  follows  immediately  that  both  vectors 
are  isotropic  lines;  that  is, 

(?■  .  (?■ . 0 
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As  the  miniinal  surface  should  real  one  has  only  to  take  Q  as  the 
complex  conjugate  of  P  and  as  the  complex  conjugate  vector 
function  of  t]. 

Exact  solution  of  simplified  differential  equation  satisfying 
shock  conditions,-  Now  since 


^y/pQ^o^  ““^x/^o^o^  -  ^1  -  (g^/ pQ^^o^) 


and  =  p^U  on  the  shock  wave^  it  follows  that 


K  X  d|^  = 


1  -  PiUi|ry/po2c^2  .F.  +  pjys^^jp^2c^2  p^u 


dy 


=  Tc(y)  dy 


(33) 


where  and  ty  are  assigned  values  according  to  equation  (l7). 

Let  the  components  of  this  vector  he  denoted  hy  [”  i(y)^ 


If  F(piUy)  is  an  analytic  function  of  y  and  possesses  an  ana¬ 
lytic  extension  in  a  domain  in  the  cr-plane^  where  cr  =  y  +  ioo^  then 
equation  (30)  is  the  differential  form 


Or, 


(34) 


dx  =  Re 


i"  (p^up) 


•f  1 


1  -  P2^Ufyjpo2c^2 


\ji  - 


dP 


(35) 
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dy  =  Re 


1  +  i 


-f  (PiW)  *  PiWjj/Po^c, 


\1  -  (b^/po^Pq^) 


dP 


(36) 


d0  =  Re 


V  (fi®)  - 


-  (b^/po^Po^) 


iPlU 


dP 


(37) 


From  the  definition  of  h(n),  hy  eqxiation  (23).  and  relation  (28),  it 
is  evident  that  the  last  eqiiation  in  equations  (27)  can  also  be  written 
as 


d0  +  i  dilf 
Po^o 


=  itf”'(t)  dt 


(38) 


Because  of  the  fact  that  ti(p)  in  equation  (31)  is  an  isotropic  line 
when  P  =  Constant,  or  (ll')2(P)  +  (n2')2(P)  .  (r|3')2(P)  =  0^  one  can 

find  an  analytic  relation  which  transforms  t  into  P  so  that 


tIj^Mp)  dP  =  i  f^"'(t)  dt  (39) 

Ti2*(P)  dP  =  dt  (40) 

and 

Ti3'(P)  dP  =  itf3_"'(t)  dt  (4l) 


18 
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Therefore  equation  (37)  can  also  be  written  as  the  following: 


d0  +  i  dilr 


The  final  solution  (26)  is  now  obtained  by  substituting  shock  condi¬ 
tions  (17)  into  equations  (35)>  (36),  and  (42). 


Massachusetts  Institute  of  Technology 
Cambridge,  Mass. ,  July  6,  1949 
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